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Resumo 



A seguinte monografia tern como objetivo o desenvolvimento e documentacao 
de urn aplicativo que auxilie no aprendizado sobre dinamica nao-linear e caos. 0 aplica- 
tivo tern como objeto de estudo o modelo do oscilador nao-linear de Duffing. Seu corn- 
portamento caotico proporciona a implementacao de diversas ferramentas interessantes 
ao estudo da nao-linearidade tais como a simulacao grafica de seu comportamento ffsico 
dentre outras. 0 aplicativo utiliza-se de metodos numericos para relacionar-se ao modelo 
ffsico do Oscilador de Duffing. Para o desenvolvimento deste trabalho foi necessario o 
estudo teorico no que concernem as diversas disciplinas envolvidas (Ffsica, Matematica, 
Estrutura de Dados, Engenharia de Software, dentre outras) e configurar os recursos e 
aplicacoes que o software, anteriormente em projeto, possui. 

Palavras-chave: Caos e Dinamica nao-linear, Oscilador de Duffing, Double- 
well, Simulador. 



Abstract 



The present monograph aims to develop and substantiate a software that helps 
on the learning of chaos and non-linear dynamics. The application has model called 
Duffing's non-linear oscillator as object of study. Its chaotic behavior allows the implement 
of several tools that are interesting for the study of non-linearity such as the graphical 
simulation of its physical behavior and others. The application rely on the use of numerical 
methods to link itself to the physical model of Duffing's oscillator. To accomplish this 
paper it was required the theorical study that envolves the several linked matters (such as 
Physics, Mathematics, Data Structure, Software Engineering and others) and the setting 
of features and applications that the software, previously in project, has. 

Key-words: Non-linear dynamics and Chaos, Duffing's oscillator, Double-well, 

simulator. 
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1 INTRODUQAO 



1.1 Motivacao 

1.2 Justificativa 

A disciplina de Oficina de Integracao I do curso de Engenharia de Computacao 
ofertado na UTFPR tern como ementa o seguinte texto: 

"Integracao dos conhecimentos de disciplinas de formacao basica e profission- 
alizante obtidos ate o momento. Aplicacao dos conceitos de metodologia cientffica para 
o desenvolvimento, em equipes, de urn sistema computational (software) contemplando 
essa integracao. Aplicacao dos conceitos de metodologia cientffica e comunicacao e ex- 
pressao para a elaboracao e apresentacao oral de relatorio final dos resultados do projeto 
desenvolvido." (UTFPR, 2006) 

Ao analisar a ementa pode-se observar que Oficina de Integracao I tern como 
objetivos a multidisciplinariedade e a integracao de outras disciplinas. Essa visao e 
reforcada no artigo "Interdisciplinariedade e Criticidade na Formacao em Engenharia 
de Computacao na UTFPR" pois, segundo os autores, "...as disciplinas de oficinas de 
integracao sao, ao mesmo tempo, estrategia e ferramenta para promover a integracao 
curricular" .(MERKLE; BEZERRA; GIMENEZ-LUGO, 2008). 

Tendo em vista este contexto, o trabalho desenvolvido se encaixa nos requisitos 
da disciplina. Pois o desenvolvimento de um aplicativo que simule o modelo do oscilador 
nao-linear de Duffing exige a integracao de varias disciplinas tais como: Ffsica, Dinamica 
nao-linear, Calculo Numerico, Fundamentos de Programacao e Engenharia de Software. 

Alem de atender os objetivos da disciplina, outra grande motivacao para a 
execucao deste projeto foi a falta de um aplicativo que retina todas as ferramentas que 
o software desenvolvido neste trabalho possui e que sao necessarias para um bom en- 
tendimento tanto de dinamica nao-linear caotica quanto do proprio oscilador de Duffing 
*****COLOCAR ALGUMA REFERENCIA DA ANALISE DE SIMILARES*****. E 
importante a diversificacao das ferramentas de ensino de Ffsica e Matematica. A mod- 
elagem computational proporciona a interacao entre os recursos audio-visuais de forma 
simples e dinamica, portanto, favorecendo o aprendizado. *********** Alguns aplicativos 
**REFERENCIA A ANALISE** foram desenvolvidos com o mesmo intuito, porem nao 
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possuem as ferramentas julgadas importantes no bom entendimento de dinamica nao- 
linear e caos. Isto fica evidente em urn relato de experiencias realizadas em sala de aula 
com o Modellus, urn software de modelagem interativa com matematica, que nao apresen- 
tou forte dependencia das condicoes iniciais do oscilador de Duffing. (MARIO BOLDO; 
CIANI, 2007). 
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1.3 Objetivos 

1.3.1 Gerais 

Auxiliar o ensino de dinamica nao-linear e caos, no estudo de modelos nao- 
lineares biestaveis caoticos, por meio de uma simulacao computacional que o software 
desenvolvido proporciona. Tornando, assim, viaveis as simulacoes de experimentos cus- 
tosos. 

1.3.2 Especificos 

Implementar algumas funcionalidades especfficas que, de acordo com a visao 
deste projeto, sao importantes para um bom aprendizado sobre dinamica nao-linear. Sao 
alguns exemplos dessas sub-rotinas: retrato ou piano de fase, secoes temporais e uma 
animacao que simule o comportamento ffsico do oscilador de Duffing que e chamada de 
espaco de configuracao. 

*** Revisao Bibliografica sobre DNL, estudo de padroes de projeto, metodos 
numericos ensino de fisica, etc. *** 



1.4 Metodologia 

1.5 Estrutura do Trabalho 
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2 ENSINO 

0 aprendizado e um processo ativo de construcao ao inves da simples recepcao 
ou aquisicao de conhecimentos (FONTOURA,2002). 

Com a difusao da computacao, cada vez mais a educacao tern como uma de 
suas ferramentas o computador. Porem, muitas vezes seu uso nao e o mais adequado. 
Quando ha computadores nas salas, muitas vezes eles sao usados da mesma maneira 
que os outros objetos mais antigos como o caderno e o quadro. Alunos usam ?Personal 
Computers? somente para fazer anotacoes enquanto professores os usam apenas como 
usariam quadros negros digitals. Torres comenta o potencial que pode ter esse tipo de 
iniciativa: 

7...0 aspecto das ajudas tecnico informaticas e enfatizado, tendo em vista que 
ele pode contribuir para uma acentuada melhoria, tanto no trabalho dos docentes como 
na aprendizagem dos alunos? (TORRES apud VIEIRA, 2002, p. 166). 

Um dos pioneiros em usar ambientes virtuais para o ensino foi Seymour Pa- 
pert, um educador matematico, que nos anos sessenta no MIT - Massachusetts Institute 
of Technology, de Cambridge, MA, Estados Unidos iniciou o projeto Logo. Adaptada 
para o portugues em 1982, na Unicamp, pelo Nucleo de Informatica Aplicada a Educacao 
(NIED), a linguagem ?Logo? vem sendo utilizada para trabalhar com criancas e adoles- 
centes. ?Logo? e uma linguagem de programacao, isto e, um meio de comunicacao entre 
o computador e a pessoa que ira usa-lo. A principal diferenca entre ?Logo? e outras lin- 
guagens de programacao esta no fato de que foi desenvolvida para ser usada por criancas 
e para que as criancas possam, com ela, aprender outras coisas. 

Uma iniciativa interessante esta no projeto Visual GD, tema de estudo da tese 
de doutorado de Josiane Wanderlinde Vieira, programa que tern como foco o ensino de 
geometria descritiva para alunos surdos. Atraves da hiperrmdia usou-se de videos, grahcos, 
texto e interatividade para ensinar um conteudo relativamente abstrato para dehcientes 
auditivos. **FAZER REFERENCIA AO PROJETO VISUAL GD** 

Um dos grandes impecilios para o aumento do uso de aplicativos computa- 
cionais nas escolas e o conservadorismo, como ahrma o professor da Unicamp Eduardo 
Chaves ?A tecnologia vem sendo introduzida na escola de forma conservadora, domesti- 
cada, o que apenas reforca o paradigma educacional tradicional, centrado na transmissao 
de informacoes disciplinares atraves do ensino do professor? (CHAVES, 2002). Porem essa 
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visao pode estar mudando. Pois ao reconhecer o computador como uma ferramenta inter- 
ativa ?...representa criar mensagens para urn novo paradigma de aprendizagem, no qual 
o aluno ao inves de ouvir e assimilar, interage ativamente com a tecnologia.? (ROMIS- 
ZOWSKI,2002) Isso entra em consenso com a afirmacao de Viet, 2002 "...a modelagem 
computacional e uma ferramenta indispensavel para o ensino de Ffsica nos dias atuais". 
Uma afirmacao extremamente importante para este estudo, uma vez que o conteudo do 
aplicativo a ser desenvolvido esta inserido no campo da Ffsica. 

2.1 Ensino de Fisica 

Ensino de Ffsica tern se alterado significativamente nas ultimas decadas devido 
a varios recursos computacionais que implementam novas ferramentas para o aprendizado 
(GLEISER, 2000). 

A modelagem e de fundamental importancia para o desenvolvimento do con- 
hecimento educacional cientffico. Esta pode ser utilizada para melhorar o entendimento 
da Ffsica. Uma das grandes crfticas das estrategias didaticas e a desconexao que es- 
tas tern entre a pratica e a teoria*** *****, de forma que a grande motivacao no estudo 
da Ffsica reside na relacao existente entre situacoes cotidianas e os modelos conceitu- 
ais. (BRANDAO; ARAUJO; VEIT, 2008) O metodo cientffico consiste na elaboracao de 
teorias que expliquem determinados comportamentos de fenomenos. A validade dessas 
teorias e veriffcada empiricamente. Todavia, na criacao de modelos cientfficos ha limi- 
tacoes e restricoes referentes a aplicabilidade desses modelos. No ensino da Ffsica, ao 
estabelecer certas idealizacoes e aproximacoes, podem-se relacionar situacoes do dia-dia 
com modelos ffsicos. Nesse contexto, a modelagem computacional permite a interacao 
entre a abordagem teorica e a pratica. (BRANDAO; ARAUJO; VEIT, 2008) 

?0 computador, visto neste processo como uma ferramenta didatica no auxflio 
da aprendizagem, pode fornecer oportunidades fmpares para a contextualizacao, visualiza- 
cao e apresentacoes das mais diversas situacoes ffsicas que possam dar sentido ao conceito 
ffsico que esteja sendo trabalhado pelo professor.? (BRANDAO; ARAUJO; VEIT, 2008) 

2.1.1 Simulacoes Computacionais no Ensino de Fisica 

O ensino da Ffsica envolve a abordagem de alguns conceitos abstratos, dos 
quais utiliza da Matematica para obter solucoes. Alem disso, a Ffsica lida com objetos de 
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estudo fora do mvel de percepcao humana, como: partfculas subatomicas, corpos dota- 
dos de grande velocidade e estudos complexos. Isto inviabiliza a representacao de certos 
modelos fisicos, pois os recursos em sala-de-aula sao limitados. Outra dificuldade esta na 
elaboracao de representacoes estaticas de fenomenos dinamicos. Alguns livros em busca 
de representar estes tipos de movimento utilizam ate fotografias estroboscopicas. As ani- 
macoes computacionais resolvem o problema de representar tais modelos fisicos pois um 
movimento, por exemplo, pode ser visualizado com extrema facilidade. As simulacoes 
computacionais oferecem maior interatividade entre os usuarios e o computador. Elas 
possuem mecanismos que permitem ao usuario inserir valores de entrada, que sao proces- 
sados, e originam animacoes graficas e outros recursos, como a elaboracao de graficos do 
movimento. (MEDEIROS,2002) 

Os modelos computacionais sao baseados em modelos fisicos reais, que sao ex- 
pressos matematicamente, processados e reconstrmdos atraves de animacoes. 0 principal 
uso da simulacao computacional e quando o experimento real e inviavel, seja por custos 
ou pela rapidez da realizacao, ou pela periculosidade do experimento ou pela interferencia 
de fatores externos. (MEDEIROS,2002) 

Sao apontados alguns beneffcios da simulacao computacional, como: 

1. Reduzir o ?rufdo? do experimento para concentrar-se somente no fenomeno ffsico; 

2. Fornecer um feedback para aperfeicoar a compreensao dos conceitos fisicos; 

3. Permitirem aos estudantes coletarem grandes quantidades de dados em menor tempo; 

4. Melhorar a compreensao dos fenomenos fisicos; 

5. Auxiliar os estudantes a aprenderem sobre o mundo natural, vendo e interagindo 
com os modelos cientfficos subjacentes que nao poderiam ser inferidos atraves da 
observacao direta; 

6. Estabelecer as relacoes de causa e efeito de fenomenos complexos; 

7. Aplicacao de conceitos abstratos; 

Entretanto, devem-se destacar algumas limitacoes impostas pelas simulacoes 
computacionais no ensino da Ffsica. Como por exemplo, que um modelo real geralmente 
e complexo, e a representacao deste em termos computacionais envolvem simplificacoes 
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e aproximacoes. Os recursos graficos das animacoes podem distanciar do modelo real e 
dificultar a compreensao de sua totalidade. As distorcoes ocasionadas da matematizacao 
do modelo fisico, se nao forem explicitadas, resultam em equfvocos fisicos. E importante 
destacar que como o modelo computacional e gerado de um modelo fisico, ambos pos- 
suem um campo de validade. Portanto, na implementacao de um software de simulacao 
fisica, e necessario considerar as restricoes e comportamentos do modelo fisico. Equfvocos 
na elaboracao de softwares de simulacoes podem comprometer o entendimento da Fisica. 
Alguns criticos afirmam que o processo de aprendizagem por erros empfricos e algo que 
nao pode ser ?simulado?. (MEDEIROS,2002) A utilizacao da ********* Informatica no 
Ensino envolve a analise critica estabelecida pelas limitacoes que a modelagem computa- 
cional favorece no ensino e suas possibilidades. ?Nao existe duvida de que a aprendizagem 
na interacao com seres humanos, com um professor real, e uma necessidade indispensavel 
em muitos casos. Simulacoes podem ser ricas em imagens e em calculos; mas sao pobres 
em sua dimensao social.? (MEDEIROS,2002) conclui-se, entao, que a simulacao com- 
putacional verificada na realizacao deste trabalho nao dispensa que seja ministrada por 
um professor em uma disciplina especffica o que sera comentado com mais detalhes na 
secao que comenta sobre "A utilizacao do software na disciplina de Dinamica nao-linear e 
Caos". 

2.2 Ensino de Dinamica Nao-Linear e Caos 

2.2.1 A utilizacao do software na disciplina de Dinamica nao-linear e Caos 

0 software simulador do Oscilador de Duffing em questao foi desenvolvido para 
posterior utilizacao na disciplina ?T6picos de Dinamica Nao-Linear e Caos? ofertada pelo 
***listas de siglas?*** CPGEI (Programa de Pos-Graduacao em Engenharia Eletrica e 
Informatica Industrial), da UTFPR, destinada aos estudantes de pos-graduacao e exten- 
siva aos estudantes de graduacao. Essa disciplina, ministrada pelo professor Mario Sergio 
Teixeira de Freitas, do DAFIS ***lista de siglas?***, tern metodologia diferenciada pois 
sua estruturacao baseia-se em aulas expositivas teoricas, com exemplos de referenda e an- 
imacoes e, em aulas praticas, nas quais os alunos desenvolvem simulacoes computacionais 
de sistemas nao-lineares e incentivado o uso de softwares interativos que demonstrem 
comportamento nao-linear. A disciplina apos ser implementada obteve elevado indice de 
aprovacao, ampla aceitacao pelos alunos e gerou trabalhos cientfficos. (FREITAS M. S. 
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T. 2008) 

3 Dinamica Nao-Linear e Caos 

Caos e o termo usado para descrever a aparente complexidade do que sao con- 
siderados sistemas simples. 0 comportamento caotico parece aleatorio, pois e influenciado 
por fatores externos e por graus de liberdade que influenciam o sistema (HILBORN,2004). 

Os sistemas naturais estudados ****** S ao determinfsticos, pois conhecendo as 
condicoes iniciais do sistema, em princfpio para dado intervalo, pode-se determinar pre- 
cisamente as suas condicoes finais. Mais precisamente, um sistema determinfstico e aquele 
que possui uma regra especffica que define o estado do sistema, ou seja, uma funcao que 
descreve os pontos do sistema. Porem, conceitualmente, parece-nos conflitante a ideia de 
um sistema de comportamento aparentemente aleatorio e, ao mesmo tempo, determinfs- 
tico. 0 que explica esse paradoxo e o conceito de nao-linearidade. Se causarmos pequenas 
perturbacoes em um sistema nao-linear, observam-se grandes divergencias na trajetoria 
de fase (KAPLAN, 1995). 0 estudo dos sistemas caoticos foi possfvel somente ha poucas 
decadas, pois exige o processamento computacional para a integracao de equacoes difer- 
enciais nao-lineares (FREITAS, 2003). Desde 1984, foram incorporados ao vocabulario 
cientffico termos como: caos, dinamica nao-linear e fractals. E importantes conceitos da 
dinamica nao-linear foram desenvolvidos e aplicados em grandes areas de estudo como: 
biologia, engenharia, medicina, ecologia, economia e astronomia (KAPLAN, 1995). 

0 caos pode ser conceituado como um comportamento dinamico aperiodico em 
um sistema determinfstico cujas alteracoes nas condicoes iniciais provocam sensibilidade 
nos resultados (HILBORN,2004). 

A propriedade aperiodica caracteriza-se por nao haver obrigatoriedade na repeticao 
das transicoes de estado de um sistema. As simulacoes computacionais dinamicas e aper- 
iodicas sao evidencias de caos. 

A sensibilidade as condicoes iniciais e uma caracterfstica fundamental para 
a existencia de caos, o que implica que a evolucao do sistema pode ser alterada por 
pequenas perturbacoes. Isso proporciona maior dihculdade em predizer resultados ao 
longo do tempo em um sistema caotico, pois as condicoes os alteram significativamente. 

Essas propriedades fazem com que o comportamento caotico possua grande 
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flexibilidade, proporciona aos sistemas intervencoes rapidas associadas a determinadas 
perturbacoes. Isto e caractenstico de sistemas naturais onde o caos e diversos mecanismos 
regulatorios fornecem uma flexibilidade incomum (SAVI,2006). 

Neste contexto, a descricao ou analise de fenomenos naturais atraves de mod- 
elos ou tecnicas nao-lineares e mais efetiva ***** do que os modelos ou tecnicas lineares. 
Contudo, as dificuldades inerentes ao estudo dos problemas nao-lineares, e o sucesso da 
mecanica linear em diversas areas do conhecimento, incentivaram cada vez mais o estudo 
de modelos lineares e bem comportados. 

Ao longo da historia, a analise nao-linear foi tradicionalmente evitada, criando 
um paradigma linear. Na citacao de Laplace, podemos verificar o determinismo estrito 
presente no paradigma linear: 

"Devemos ver o estado presente do universo como o efeito do seu estado ante- 
rior, e como a causa daquele que vira. Uma inteligencia que, em qualquer instante dado, 
soubesse todas as forcas pela qual o mundo natural se move e a posicao de cada uma 
de suas partes componentes, e que tivesse tambem a capacidade de submeter todos estes 
dados a analise matematica, poderia encompassar na mesma formula os movimentos dos 
maiores objetos do universo e aqueles dos menores atomos; nada seria incerto para ele, e 
o futuro, assim como o passado, estaria presente diante de seus olhos." (SAVI,2006). 

No final do sec. XIX Henri Poincare, matematico e ffsico trances, comecou a 
estudar a estabilidade do Universo. Este estudo viria a se contrapor com o determinismo 
estrito de Laplace, como pode-se ver na citacao de Poincare: ?Uma causa muito diminuta, 
que nos escapa, determina um efeito consideravel, que nao podemos deixar de ver, e entao 
dizemos que este efeito e devido ao acaso. Se pudessemos conhecer exatamente as leis da 
natureza e a situacao deste mesmo universo no instante subsequente, mesmo quando as 
leis naturais ja nao tivessem mais segredo para nos, so poderfamos conhecer a situacao 
inicial aproximadamente. Se isto nos permite aproximar a situacao subsequente com o 
mesmo grau de aproximacao, ficamos satisfeitos, dizemos que o fenomeno foi previsto, 
que e governado por leis. Mas nem sempre isto ocorre; pode acontecer que diferencas 
mfnimas nas condicoes iniciais produzam diferencas muito grandes no fenomeno final; um 
erro mmimo nas primeiras produziria um erro enorme nesses ultimos. A previsao torna-se 
impossfvel e temos o fenomeno ao acaso.? (SAVI,2006) 

Apesar de Poincare ter uma visao clara sobre a sensibilidade do sistema as 
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condicoes iniciais, foi so em 1963 com Lorenz, que houve uma analise mais aprofundada 
de sistemas dinamicos. Ele desenvolveu estudos sobre problemas atmosfericos. Contando 
com o auxilio de um computador, Lorenz tratava o modelo de Rayleigh-Benard para a 
conveccao de fluidos, e observou que uma pequena variacao nas condicoes iniciais pode- 
ria acarretar grandes diferencas na evolucao do sistema (SAVI,2006). 0 comportamento 
caotico nao-linear do oscilador de Duffing pode ser observado com os retratos de fase do 
movimento que sao graficos em que, no eixo das ordenadas, se encontram os valores de 
velocidade e nas abcissas, valores de posicao. Os parametros para a analise desse movi- 
mento sao: o coeficiente de amortecimento viscoso do meio, a constante de elasticidade da 
haste metalica, a atracao magnetica dos imas, a frequencia e a amplitude da oscilacao de 
fonte externa****REFERENCIA AO HOLMES SOBRE OS PARAMETROS****. Ini- 
cialmente, o usuario fornece os dados iniciais da velocidade e da posicao da partfcula. 
0 sistema modelado matematicamente com a utilizacao de equacoes diferenciais. E sua 
resolucao envolve metodos de integracao numericos. 0 estudo do sistema de Duffing nao 
e de grande complexidade, pois envolve osciladores que podem ser tratados unidimension- 
almente e as estruturas fractals sao identfficadas com a ?secao de Poincare?. A ?secao 
de Poincare? e um modo de detectar solucoes periodicas para a observacao de um corte 
transversal (FREITAS, 2003). Denomina-se tambem ?amostragem estroboscopica? nos 
sistemas dependentes do tempo. A simulacao numerica e a interpretacao geometrica do 
piano de fase sao recursos que podem ser utilizadas desde a sistemas simples, como o os- 
cilador harmonico simples, a sistemas complexos, como osciladores com excitacao externa 
e amortecimento viscoso. (FREITAS, 2003) 



3.1 Oscilador de Duffing 
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Figura 1: 0 Oscilador de Duffing 
A equacao diferencial que caracteriza o oscilador de Duffing pode ser descrita 

como: 

3^ x doc 

— - = -x 3 + ax - 2(5 h Acos(ojt) (1) 

dt 2 dt 

A Figura ?? ilustra o Oscilador de Duffing, que consiste em uma haste metalica de 
elasticidade a fixada a um suporte que oscila com amplitude A e freqiiencia angular uo, 
tal oscilacao realiza o papel de uma forca excitadora externa. 0 oscilador esta sujeito, 
alem da forca mensionada, ao amortecimento viscoso f3 do meio em que esta inserido e 
a acao de dois mias identicos eqiiidistantes da posicao da haste. F ilustrado, tambem, o 
ponto estudado que se move segundo um arco de cirfunferencia, porem e adotado que este 
se move segundo um eixo x em uma direcao horizontal. Esta aproximacao pode ser feita 
pois a amplitude de oscilacao deste ponto e pequena comparada ao tamanho da haste. 
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A funcao da energia potencial deste oscilador e descrita como: 

U(x) = ~kx 2 + jj-x 4 (2) 



Baseado em experiencias com o oscilador harmonico simples, sabe-se que a 
partfcula estudada reagira com um movimento de amplitude relativamente grande quando 
a frequencia da forca externa coincide com a frequencia de oscilacao natural da mesma 
partfcula. 0 problema e que a frequencia de oscilacao natural depende da amplitude 
deste movimento. Entao como a partfcula comeca a reagir a forca excitadora externa, 
a amplitude muda e por isso, a frequencia de oscilacao natural tambem muda. Varias 
caracterfsticas singulares podem aparecer (HILBORN, 2004, p. 581): 

1. A curva de resposta do sistema muda com o aumento da amplitude da forca excita- 
dora externa; 

2. The response curve of the system shows hysteresis: the response amplitude depends 
on whether we increase the frequency through the resonance region or decrease the 
frequency through the resonance region. 

3. 0 sistema pode apresentar comportamento caotico. 



3.1.1 Graficos 

• Piano de Fase: 0 piano de fase e um grahco onde, no eixo das abcissas, constam os 
valores da posicao da ponta da haste do oscilador e, o eixo das ordenadas, constam 
os valores da velocidade deste mesmo ponto. 

• Secoes Temporais: As secoes temporais sao grahcos onde o eixo das abcissas cor- 
responde ao tempo e o eixo das ordenadas e o dado que se pretende analizar por 
exemplo: Grahco da posicao contra o tempo, velocidade contra o tempo e da acel- 
eracao contra o tempo. 
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• Secao de Poincare: A secao de Poincare, e obtida atraves de uma analise estro- 
boscopica no piano de fase, ou seja, registramos a posicao da ponta da haste do 
oscilador em urn dado instante inicial e repitimos isso periodicamente em um certo 
intervalo de tempo adotado. Veremos que as posicoes nesses instantes formam o que 
e chamado de Secao de Poincare. 

• Grafico de Energia: 0 esboco do grafico da energia, pode ser obtido pela plotagem 
da equacao (2) e o grafico desta funcao e chamado de double-well. 0 nome ?double- 
well? esta relacionado a correspondencia entre a funcao da energia potencial com a 
estrutura da mola (HILBORN,2004). 



Figura 2: Piano de Fase 



Figura 3: Segao de Poincare 
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Figura 4: Grdfico de Energia 



4 Software 

Segundo PRESSMAN, 1995 software pode ser caracterizado como: ?(1) in- 
strucoes(programas de computador) que, quando executados, produzem a funcao e o 
desempenho desejados; (2) estruturas de dados que possibilitam que os programas ma- 
nipulem adequadamente a informacao; e (3) documentos que descrevem a operacao e o 
uso dos programas?. Ainda citando o mesmo autor, ele tambem define tres caracterfsticas 
que distinguem softwares de hardwares: 

1. 0 software e desenvolvido ou projetado por engenharia, nao manufaturado no sen- 
tido classico; 

2. Software nao se ?desgasta?; 

3. A maioria dos softwares e feita sob medida em vez de ser montada a partir de 
componentes existentes. 

Pressman divide os softwares em sete grupos diferentes: Basico, de Tempo 
Real, Comercial, Cientffico e de Engenharia, Embutido, de Computacao Pessoal, de In- 
teligencia Artificial. 

0 software a ser desenvolvido neste projeto encaixa-se no grupo Cientffico e 
de Engenharia. Este grupo tern sido caracterizado por algoritmos de processamento de 
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numeros. Suas aplicacoes variam desde astronomia ate vulcanologia. 

4.0.2 Engenharia de Software 

Engenharia de Software e o estabelecimento e uso de solidos princfpios de 
engenharia para que se possa obter economicamente um software que seja confiavel e 
que funcione eficientemente em maquinas reais (PRESSMAN, 1995). Alem do aspecto 
economico esta disciplina atraves de seus metodos e ferramentas possibilita um desen- 
volvimento racional e oportuno do software de computador (PRESSMAN, 1995). 
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4.1 Analise de Aplicativos Similares 

A Figura 7 apresenta uma analise que tern como objetivo a busca e o levanta- 
mento dos pontos fortes e fracos de aplicativos similares ao proposto neste projeto. Este 
tipo de estudo permite que o projeto tenha tanto uma base para seu imcio quanto uma 
base de comparacao. 0 estudo foi tabulado, sendo as colunas representando os aplicativos 
e as linhas os elementos avaliados. A variavel de avaliacao tern tres estados possfveis como 
se pode ver na Figura 6. 



Figura 5: Legenda da Analise 

A classificacao positivo e dada quando a ferramenta avaliada existe no aplica- 
tivo e sua configuracao e satisfatoria. Ja negativo e atribufda quando a ferramenta existe 
mas sua configuracao e confusa ao ponto de poder ser considerada nao-satisfatoria. E por 
ultimo a classificacao nao possui quando o aplicativo nao disponibiliza esta ferramenta. 
A seguir os enderecos eletronicos onde foram encontrados os aplicativos em questao e a 
tabela comparativa. 

• A - http://www.math.udel.edu/~hsiao/m302/JavaTools/osduffng.html 

• B - http://www.maths.usyd.edu.aU/u/nigel/NLD/Duffing/Duffing.html 

• C - http : //www . compadre . org/OSP/document/ServeFile . cf m?ID=7857&DocID= 
642 

• D - http://www.df.ufcg.edu.br/~adriano/ClickApp2.html 

• E - http : //brain. cc . kogakuin . ac . jp/~kanamaru/Chaos/e/Duf f ing2/ 

• F - Software desenvolvido neste trabalho. 



gura 6: Andlise 
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4.2 Padroes de Projeto 

Este trabalho utiliza o que e chamado de Padroes de Projeto ou Design Pat- 
terns. Estes padroes sao ferramentas que permitem um melhor desenvolvimento do aplica- 
tivo, pois eles solucionam muitos dos problemas que os projetistas enfrentam no dia a dia 
(Gamma et al,2000 p.27).Outro beneffcio que estas ferramentas trazem esta na possibili- 
dade de reutilizacao de solucoes, algo que pode acelerar o desenvolvimento do aplicativo. 

4.2.1 Padrao MVC 

A sigla MVC vem do ingles (Model, View, Controller) que pode ser traduzido 
como Modelo, Visao e Controle. Segundo Lopes, 2006 MVC e um padrao de arquitetura 
de Software criado na linguagem Smalltalk e tern como objetivo desacoplar a interface 
grafica do codigo que realizara o trabalho. Abaixo uma breve descricao de como atua 
cada parte desta arquitetura: 

• Visao: tern como funcao mostrar dados e receber informacoes adicionadas pelo 
usuario. Controle: faz a comunicacao entre os dados recebidos pela interface e 
o modelo. Tambem nesta camada ocorre o processamento da informacao que sera 
enviada ao modelo. 

• Modelo: e a parte do software que ira processar dados e que executara as funcoes, 
calculos, metodos necessarios para cumprir o objetivo do programa. 

Um dos padroes de projetos utilizados neste tipo de abordagem e o padrao observer 
(observador) que e definido por Gamma como o padrao que ?.. define uma dependencia 
um-para-muitos entre objetos para que quando um objeto mudar de estado, todos os 
seus dependentes sejam notificados e atualizados automaticamente."Este padrao e muito 
importante neste projeto, uma vez que varias das ferramentas do software irao mostrar 
informacoes provenientes de uma mesma fonte, mas de formas diferentes. Alem disso 
o uso desse padrao permite uma expansao do sistema de forma sustentavel, pois se for 
preciso adicionar uma nova ferramenta, nao sera necessario mexer no cerne do programa. 
(FERNANDES,2006) 
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4.2.2 Importancia do Padrao de Projeto Observer no Sistema 

0 padrao observer tern um grande destaque neste trabalho, pois o seu fun- 
cionamento e uma das pecas cruciais para que o sistema possa propiciar a experiencia 
objetivada ao usuario. Isto ocorre pois este padrao permite que varios objetos sejam 
atualizados automaticamente de acordo com mudancas no modelo. Sua principal apli- 
cacao esta na parte de visualizacao de dados em que sao necessarias varias representacoes 
diferentes dos mesmos, o que e o caso deste trabalho. Este padrao tern como dois atores 
uma interface chamada Observador e uma classe abstrata Observado. A classe Obser- 
vado conhece todas as classes que implementam a interface Observador. Assim quando 
o Observado e atualizado, ele notifica todos os observadores de sua mudanca. Estes por 
sua vez recebem os novos dados e os tratam da maneira que precisam para executar sua 
funcao. No caso deste projeto todas as ferramentas implementam a interface Observador 
e o modelo deriva da classe Observado, o que significa que quando houver uma atualizacao 
no modelo, todos as ferramentas serao atualizadas. 

Exemplo de funcionamento: 



Figura 7: Exemplo de Funcionamento do Observer 
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0 metodo setDados da classe Modelo adiciona o novo tempo a funcao de 
integracao fazendo com que o modelo tenha os proximos valores do sistema. Alem disso 
ele chama o metodo notificar herdado da super classe abstrata Observado. Esse metodo 
por sua vez chama em todos os observadores o metodo update. Este metodo tern sua 
implementacao especffica em cada ferramenta juntamente para que cada um trate os 
dados da maneira conveniente. 0 codigo acima mostra o exemplo dessa comunicacao 
para a classe da JanelaPlanoDeFase e JanelaMostrador. 

4.3 UML 

A Unified Modeling Language e atualmente o esquema de representacao gra- 
fica mais amplamente utilizado para modelar sistemas orientados a objetos.(DEITEL e 
DEITEL, 2005). 

4.4 JAVA 

A linguagem Java foi escolhida para o desenvolvimento deste projeto por ser 
uma linguagem geralmente utilizada com o paradigma de orientacao a objetos, alem de 
sua grande popularidade. Esta linguagem tern mais de 6.5 milhoes de desenvolvedores . A 
popularidade da linguagem assegura um grupo imenso de recursos de suporte (BARNES; 
ROLLING, 2004). E quanto mais recursos mais faceis sera superar os obstaculos que 
podem vir a ser encontrados durante a fase de desenvolvimento do aplicativo. Alem destes 
motivos tambem foi levado em conta que e a linguagem de programacao comum a todos 
os integrantes do projeto, assim aproveitando ao maximo os conhecimentos ja adquiridos 
pelos componentes da equipe. Outro ponto a favor desta linguagem e sua vasta API 
(aplication programmng interface), que possui varios conjuntos de classes prontas que o 
desenvolvedor pode utilizar-se. 

4.5 Interface Grafica 

0 sistema desenvolvido possui uma grande gama de ferramentas, sendo o 
usuario a decidir quais sao relevantes no momento em que esta interagindo com o pro- 
grama. Para propiciar essa liberdade, optou-se por uma abordagem com janelas internas 
ao programa, sendo elas facilmente fechadas e abertas. 0 projeto das janelas deste pro- 
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grama difere urn pouco de outros, uma vez que normalmente este tipo de janela pode ser 
fechado, minimizado e aberto. Neste caso, como as janelas sao predefinidas elas nao podem 
ser minimizadas, pois nao existe uma funcao relevante para tal tipo de comportamento. 
Assim o comportamento basico das janelas e aberto ou fechado dentro do pro grama. As 
janelas ativas ficam em uma area grande do aplicativo e existe uma outra area menor onde 
ficam os fcones que representam as janelas que nao estao em uso. Essa divisao existe para 
que o usuario tenha um acesso rapido a todas as ferramentas que desejar, pois caso clicar 
no fcone referente a ferramenta desejada, esta ira abrir para ser utilizada. A maioria das 
ferramentas sao grahcos que auxiliam no estudo do movimento caotico executado pelo 
sistema. 

4.5.1 Icones 

Os fcones das ferramentas possuem cada um uma cor, uma imagem representa- 
tiva e o nome da ferramenta. Esses tres tipos de informacao tern como objetivo diferenciar 
uma ferramenta da outra, tentando evitar confusao. 

4.6 Analise Orientada a Objetos 

4.6.1 Descricao do Problema 

Para que o software desenvolvido atenda aos objetivos iniciais deste projeto, 
deve propiciar a todos seus usuarios uma experiencia satisfatoria. Os atores envolvidos 
sao basicamente de dois tipos podendo haver desdobramentos, mas estes nao devem fugir 
muito dos perhs iniciais. Sao eles: 

• Professores Este tipo de usuario deve usar o sistema como uma ferramenta de ensino 
ou para alguma atividade similar. Ele deve ter acesso a todas as funcionalidades e 
o programa deve atender a necessidade explanativa que uma aula possa a vir exigir. 

• Alunos/Interessados Este ator difere do anterior em alguns fatores. Neste grupo 
existe a possibilidade de que o usuario nao seja um profundo conhecedor do tema, e 
essa limitacao nao deve o afastar do produto deste projeto. Outro fator importante 
e o carater explorativo que o software pode propiciar, fazendo com que para aqueles 
que nao tern conhecimento nenhum, ao menos possam se familiarizar com elementos 
do objeto de estudo. 
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4.6.2 Requisitos Funcionais do Sistema 

As funcionalidades de software sao as mesmas para os dois atores. 
Elementos para o estudo do sistema conceitual: 

• Grafico da posicao contra o tempo; 

• Grafico da velocidade contra o tempo; 

• Grafico da aceleracao contra o tempo; 

• Grafico da velocidade contra a posicao (Piano de Fase); 

• Grafico de Energia do Sistema; 

• Controles para interacao com o usuario; 

• Controle da constante elastica da barra; 

• Controle da constante de amortecimento do meio; 

• Controle da amplitude da forca excitadora; 

• Controle da frequencia angular da forca excitadora; 

• Controle da velocidade inicial do movimento; 

• Controle da posicao inicial do movimento; 

• Pausar a simulacao; 

• Entrada de dados numericamente (por input box); 
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4.6.3 Diagrama de Classes 



Figura 8: Diagrama de Classes 
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4.6.4 Diagramas de Sequencia 



Figura 9: Diagrama de Entrada de Dados 



gura 10: Diagrama de Agoes do Usudrio 
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4.6.5 Dicionario de Classes 

• Controle 

Classe que reuni os comandos que o usuario pode executar na interface grafica, e os 
repassa para o modelo. 

• Modelo 

Classe que executa a integracao da equacao diferencial deste estudo, armazenando 
e devolvendo os valores a serem mostrados pelas ferramentas do programa. 

• Observado 

Classe abstrata com metodos de comunicacao com objetos que implementem a in- 
terface observador. 

• Observador 

Interface que cria a relacao com a classe abstrata Observado. Essa relacao tern com 
intuito a atualizacao das ferramentas sendo usadas no sistema. 

• PrincipalJanela 

Classe que cria a janela principal do aplicativo. Derivada da classe JFrame. 

• DesktopPrincipal 

Classe responsavel pela area em que as ferramentas ativas sao mostradas. Derivada 
da classe JDesktopPane 

• Painellnferior 

Classe responsavel pela area de janelas desativadas. Derivada da classe JPanel. 

• Botaolnter 

Classe que cria o botao que representa as janelas desativadas. Derivada da classe 
JButton 

• Janelalnterna 

Super classe de to das as janelas internas. Nela estao os comportamentos comuns a 
todas as ferramentas. Derivada da classe JinternalFrame e implementa a interface 
Observador 
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JanelaConfig 

Classe que mostra a conhguracao ffsica do objeto simulado. Derivada da classe 
Janelalnterna. 

JanelaEnergia 

Classe que representa a energia do sistema. Derivada da classe Janelalnterna. 
PainelEnergia 

Classe que desenha o grahco da energia do sistema. Derivada da classe JPanel. 
JanelaEntradaDados 

Classe que em que pode-se inserir valores para alterar o comportamento do sistema. 
Derivada da classe Janelalnterna. 

PainelEntradaDados 

Classe que desenha botoes, campos, legendas referentes a ferramenta de entrada de 
dados. Derivada da classe JPanel. 

JanelaGrafico 

Super classe dos grahcos com uma variavel contra o tempo. Possui todo o compor- 
tamento padrao deste tipo de grafico no sistema. Derivada da classe Janelalnterna. 

PainelGrafico 

Classe que desenha o grahco da classe JanelaGrahco. Derivada da classe JPanel. 
JanelaXT 

Classe que representa a ferramenta referente ao grahco do espaco pelo tempo. 
Derivada da classe PainelGrahco. 

JanelaVT 

Classe que representa a ferramenta referente ao grahco da velocidade pelo tempo. 
Derivada da classe PainelGrahco. 

JanelaAT 

Classe que representa a ferramenta referente ao grahco do espaco pelo tempo. 
Derivada da classe PainelGrahco. 
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• JanelaMostrador 

Classe que apenas mostra os valores obtidos pelo modelo. Derivada da classe 
Janelalnterna. 

• PainelMostrador 

Classe que desenha na tela os valores referentes a classe JanelaMostrador. Derivada 
da classe JPanel. 

• JanelaPlanoDeFase 

Classe que representa a ferramenta referente ao grahco da velocidade pelo espaco. 
Derivada da classe PainelGrafico. 

• PainelPlanoDeFase 

Classe que desenha o grafico da classe JanelaPlanoDeFase. Derivada da classe 
JPanel. 

• JanelaPoincare 

Classe que representa a ferramenta referente ao grafico da velocidade pelo espaco. 
Derivada da classe PainelGrafico. 

• PainelPoincare 

Classe que desenha o grahco da classe JanelaPoicare. Derivada da classe JPanel. 

5 Consideracoes Finais 

5.1 Discussao 

Como resultado desta pesquisa temos um prototipo do aplicativo proposto 
inicialmente. Esse prototipo possui todos os requisitos funcionais levantados no comeco 
do projeto exceto o esboco do grahco da secao de Poincare efetivamente pois por alguma 
incompatibilidade de natureza numerica o corte de Poincare e plotado em meio a uma 
poeira de pontos o que dihculta a visualizacao, mas nao impossibilita que tal etapa seja 
atingida apos varios testes ahm de encontrar uma conhguracao mais adequada. Testes 
com usuarios se fazem necessarios para validar a ehciencia da abordagem escolhida. Assim 
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um proximo passo para este projeto seriam estas avaliacoes (testes informais foram feitos 
durante a producao do aplicativo, mas nao documentados) e mais importante ainda, 
testes em salas de aula. Avaliacoes em sala de aula seriam de grande interesse pois 
observacao em campo pode levantar inumeras situacao nao previstas pelos envolvidos 
no projeto, o que ajudaria ainda mais na evolucao deste software. A versao prototipo 
do Simulador Interativo do Oscilador de Duffing foi executada em versao para desktop. 
Outro possfvel aperfeicoamento seria uma versao online do software. Porem este tipo de 
adaptacao exigiria um outro tipo de conhecimento tecnico e teorico, pois a experiencia 
que o hipertexto proporciona e diferente da proporcionada por um aplicativo offline. 

5.2 Resultados 

0 sofware desenvolvido e composto por nove janelas internas (como se pode ver 
na Figura 12), que em conjunto simulam o comportamento caotico nao-linear do Oscilador 
de Duffing. A janela interna de Entrada de Dados permite a insercao de parametros adi- 
cionais no estudo desse movimento, a barra deslizante proporciona ao usuario o controle 
de velocidade de simulacao do software. 0 aplicativo possibilita a visualizacao da sensi- 
bilidade as condicoes iniciais do oscilador. 0 Espaco de Configuracao e uma animacao 
da ponta haste metalica se deslocando em uma dimensao. 0 grafico de Energia Potencial 
apresenta comportamento regular, pois a equacao de Energia Potencial esta em funcao 
unicamente da posicao e esta varia dentro de um intervalo definido de valores. As secoes 
temporais, posicao por tempo, velocidade por tempo e aceleracao por tempo apresentam 
comportamento nao-periodico dependendo dos valores iniciais e dos parametros que sao 
inseridos. 0 sofware proporciona um redimensionamento da maioria das janelas internas 
para melhor visualizacao e possui botoes de abertura e fechamento das mesmas. 0 Piano 
de Fase refere-se a plotagem da velocidade pela posicao. A forca excitadora do oscilador 
mantem frequencia constante, a cada iteracao proveniente de um periodo completo desta 
forca e plotado, em um piano de fase, a posicao e a velocidade neste instante. 0 resultado 
dessa iteracao da origem a Secao de Poincare. 
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Figura 11: Programa com as Janelas Internas Abertas 



Figura 12: Programa com as Janelas Internas Fechadas 
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5.3 Conclusao 

0 trabalho desenvolvido baseado nas propostas anteriormente apresentadas 
atingiu grande parte do objetivo proposto pois podemos simular o oscilador nao-finear bi- 
estavel caotico de Duffing e, por conseguinte, identificamos o seu comportamento caotico, 
por exemplo, atraves de sua sensibilidade as condicoes iniciais. Nesse contexto, a cri- 
acao de um simulador do Oscilador de Duffing contribui para a aplicacao de conceitos de 
nao-linearidade e caos. E relevante que o estudante de ffsica, mesmo que com recursos 
limitados, possa correlacionar aspectos teoricos com praticos (GLEISER, 2000). 
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Apendice 

1 Um breve manual de utilizacao 

1.1 Janelas Internas 

• Grafico X t 

Apresenta o grafico da secao temporal em que contem valores de posicao no eixo 
das ordenadas e valores de tempo no eixo das abcissas; 

• Grafico V t 

Apresenta o grafico da secao temporal em que contem valores de velocidade no eixo 
das ordenadas e valores de tempo no eixo das abcissas; 

• Grafico A t 

Apresenta o grafico da secao temporal em que contem valores de aceleracao no eixo 
das ordenadas e valores de tempo no eixo das abcissas; 

• Grafico E x 

Apresenta o grafico em que contem valores de energia potencial no eixo das orde- 
nadas e valores de posicao no eixo das abcissas; 

• Dados Analiticos 

Apresenta o valor instantaneo de posicao, velocidade e aceleracao; 

• Piano de Fase 

Apresenta o Piano de Fase do oscilador, ou seja, o grafico que apresenta valores de 
velocidade no eixo das ordenadas e valores de posicao no eixo das abcissas; 

• Secao de Poincare 

Apresenta, em um piano de fase, um prototipo da Secao de Poincare; 

• Entrada de Dados 

Abrira uma janela interna que ira conter: 
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— Pos Inicial 

Neste campo, o usuario pode configurar a posicao inicial da ponta da haste 
metalica do oscilador. Para resultados satisfatorios, e recomendado que o 
usuario coloque valores no intervalo de -2.0 a 2.0. 

— Vel Inicial 

0 usuario pode configurar a velocidade inicial da ponta da haste metalica do 
oscilador. Para resultados satisfatorios, e recomendado que o usuario coloque 
valores no intervalo de -2.0 a 2.0. 

— Coef Elast 

Neste campo o usuario altera o coeficiente elastico da haste metalica. 

— Amort Vise 

Neste, o usuario pode conhgurar o amortecimento viscoso do meio em que o 
oscilador esta inserido. 

— Amplitude Fexc 

0 usuario conhgura, neste campo, a amplitude de movimento que a forca har- 
monica excitadora desempenha. 

— Freq Fexc 

Neste campo, e conhgurado a freqiiencia da forca harmonica excitadora. 

— Botao Iniciar 
Inicia simulacao. 

— Botao Pausar 
Pausa simulacao. 

— Barra deslizante inferior 

Controla a velocidade de simulacao. Vale a pena ressaltar que esta veloci- 
dade depende da maquina em que o programa esta sendo executado e que a 
quantidade de janelas interiores abertas tambem influem neste desempenho. 

• Espaco de Conhguracao 

Abrira uma janela que ira conter a ponta da haste metalica representada por um 
cfrculo vermelho se movendo na direcao do eixo das abcissas. Contera, tambem, 
dois quadrados azuis que representam a posicao de cada una. 
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7.1.2 Abrindo Janelas Internas 

0 usuario pode escolher a janela interna que gostaria de visualizar. Para isto 
basta posicionar o ponteiro do mouse em algum dos botoes localizados na parte inferior 
da janela principal (Figura 13) e verificar, logo acima, o texto que este contem e, com 
isso, selecionar a janela desejado. 

7.1.3 Iniciando Simulacao 

Com o aplicativo aberto, assim como na Figura 13 o usuario podera clicar, com 
o ponteiro do mouse, encima do penultimo botao localizado na parte inferior da janela 
principal no sentido da esquerda para a direita. Assim abrira a janela interna Entrada 
de Dados onde o usuario podera iniciar a simulacao no botao Iniciar bem como pausar 
quando quiser clicando no botao Pausar. 



